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I.   INTRODUCTION 

The  representation  of  terrain  in  many  Army  combat 
simulation  models  is  both  costly  and  time  consuming.  For 
instance,  terrain  input  for  high  resolution  models  such  as 
Dynamic  Tactical  Simulation  (DYNTACS) ,  Tank  Exchange  Model 
(TXM) ,  and  various  versions  of  the  Bonder/Individual  Unit 
Action  model  comes  from  computer  tapes  prepared  by  such 
agencies  as  the  Army  Map  Service  and  the  Waterways 
Experiment  Station.  These  tapes  typically  contain  terrain 
data  at  ten  meter  intervals  for  an  entire  piece  of  actual 
terrain.  To  survey  and  record  the  18,000  grid  intersections 
for  a  six  kilometer  by  three  kilometer  section  of  terrain  is 
expensive  andj  as  is  suggested  by  this  paper,  unnecessary 
for  many  applications.  The  storage  of  this  information  not 
only  can  reguire  a  major  portion  of  a  computers  storage 
capacity  (core) ,  but  also  require  a  great  deal  of  running 
time  searching  this  stored  data  in  order  to  make  movement 
and  line-of-sight  calculations.  The  exact  amount  of  core 
and  running  time  depends  upon  the  type  computer,  the 
resolution  of  the  model,  and  the  dimensions  of  the  selected 
peice  of  terrain. 

The  current   method   of  terrain  representation  is  also 

statistically   questionable.     In    a    Vector    Research, 

Incorporated    (VRI)   study   on   terrain   line-of-sight,   it 

concluded  that  "...present   and   past   Army   study   results, 

based   on   the  analysis  of  combat  results  on  a  very  limited 

sample  of  terrains,  may  have  been  determined  by  the   terrain 

selection   process   and   not   by  the  actual  weapon  system  or 

force  design  differences. "   More  specifically  it  stated  the 

following: 

a.  "There  is  extreme  sensitivty  in  combat  model 
results  as  the  scenarios  (terrain  and  movement 
assumptions)  are  varied,  even  when  variation  is 
within  a  class  of  scenarios  chosen  for  their  "a 
priori"  equivalence. 


b.  "This  sensitivity  can  be  slightly  reduced,  but 
remains  extreme  (with  probabilities  of  win 
estimable  only  within  plus  or  minus  25%)  even  when 
battle  results  are  used  to  redesign  scenarios." 

These  results  imply  that  sufficient  replications  of 
each  type  of  terrain  should  be  run  in  order  to  reach  a 
satisfactory  statistical  level  of  significance.  For  VRl»s 
analysis  at  least  50  replications  of  each  type  of  terrain 
were  used. 

The  problems  stated  above  suggest  that  new  methodology 
for  the  representation  of  terrain  is  necessary  in  order  to 
reduce  costs,  shorten  computer  running  time,  and  improve  the 
level  of  significance  for  Army  studies  which  use  combat 
models.  Consequently,  the  objective  of  this  paper  is  to 
develop  this  new  technology  using  parameterized,  randomly 
created  terrain,  and  to  verify  the  concepts  using  a 
computerized  terrain  simulation  model. 


"H  Farrell,  Robert  L.,  Freedman,  Richard  J., 
Investigations  of  the  Variation  of  Combat  Model  Predictions 
witn  "Terrain  Line  oTT~  "5 igU t ,  p.  7T  vector  "Research, 
Incorporated,  1975. 

2.   Ibid.,  p.  6. 
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II.   TERRAIN  REPRESENTATION 

A.   STATE  OF  THE  ART 

DYNTACS,  which  is  one  of  the  most  frequently  used  high 
resolution  combat  simulations,  serves  as  a  good  example  of 
the  state  of  the  art  in  terrain  representation  and  use.  In 
general  terms  DYNTACS  is  a  two-sided,  dynamic,  Monte  Carlo, 
highly  interactive  combat  simulation  capable  of  representing 
forces  from  a  single  crew  served  weapon  to  a  reinforced 
battalion.  In  particular  the  model  considers  in  detail  the 
effects  of  terrain  on  detection,  mobility,  tactics,  and 
intervisibility  between  weapon  elements.  For  the  remainder 
of  this  paper,  unless  otherwise  specified,  only 
macro-terrain  (elevation  and  slope)  will  be  considered.  As 
previously  mentioned  this  data  is  supplied  by  the  Waterways 
Experiment  Station  (WES)  on  computer  tapes.  Although  these 
tapes  have  a  ten  meter  resolution,  DYNTACS  is  generally  run 
using  100  meter  squares.  The  simulation  divides  each  grid 
square  diagonally,  thus  producing  a  series  of  adjoining 
triangular  terrains.  The  entire  battlefield,  therefore,  is 
represented  as  a  surface  of  equally  sized  triangles  which 
vary  in  slope  depending  upon  the  elevation  at  their  corners. 

From  the  macro-terrain  data,  line-of-sight  between  any 
two  opposing  elements  is  computed.  This  is  accomplished  by 
first  computing  an  angle  between  the  horizontal  and  a 
straight  line  drawn  between  an  observer  and  target  (O-T 
line) .  The  program  then  conducts  a  search  of  the  terrain 
along  the  path  of  the  O-r  line  to  see  if  any  macro-terrain 
is  higher  than  the  0-T  line  itself.  This  is  accomplished  by 
comparing  the  angle  of  the  0-T  line  with  the  angle  above 
horizontal  cf  the  observer-terrain  line.  If  the  latter 
angle  is  larger,  there  is  no  intervisibility.  Over  the 
duration  of  a  battle  with  numerous  elements  this  calculation 
may  be  made  thousands  of  times.  Consequently,  not  only  is 
the   time   tc   prepare   the  tapes  high,  but  also  the  time  to 


compute  lines-of-sight  once  the  terrain  is  input  to  the 
model.  Considerable  effort  by  developers  and  users  of  this 
model  has  been  expended  to  streamline  this  subroutine. 

B.   A  DIFFERENT  APPROACH 

If  a  combat  simulation  is  to  be  used  for  determining  who 
will  win  a  battle  over  a  specific  piece  of  ground  such  as 
the  Fulda  Gap,  then  that  terrain  should  be  modelled  with 
precision.  However,  if  the  model  is  to  be  used  to  determine 
weapons  effectiveness,  then  representation  of  a  specific 
terrain  is  not  necessary,  and  need  only  be  representative  of 
the  area  to  be  modelled.  Furthermore,  as  VRI  suggested,  a 
single  realization  (e.g.  the  Fulda  Gap)  may  produce 
inaccurate  conclusions  in  a  weapons  system  evaluation.  What 
is  needed,  therefore,  is  random  terrain  representation  which 
can  be  replicated  quickly.  However,  at  the  same  time  it 
must  be  representative  of  a  particular  "type"  of  terrain. 
One  approach  is  to  represent  hills  precisely  using  a  series 
of  mathematical  equations.  This  method  would  eliminate  the 
inaccuracies  of  digitized  (discrete)  terrain  but  would  still 
have  the  problem  of  producing  results  based  on  a  single 
realization  cf  terrain.  Moreover,  the  time  involved  to 
model  with  precision  would  negate  its  contribution.  At  the 
other  end  of  the  terrain  modelling  spectrum  is  to  use 
cylinders  or  cones  to  represent  hills.  This  is 
mathematically  appealing,  but  lacks  realism  even  if  the 
method  were  proven  statistically  satisfactory.  The 
solution,  therefore,  lies  somewhere  in  between. 

An  equation  that  produces  what  appears  to  be  a  simple 
hill  and  at  the  same  time  is  mathematically  tractable  is  the 
bivariate  normal  probability  density  function.  Figure  1  is 
a  sample  realization  of  terrain  created  and  drawn  by  a 
computer,  using  the  bivariate  normal  distribution.  The 
second  figure  is  a  contour  map  of  the  same  terrain,  also 
drawn  by  the  computer. 
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Figure  2:   The  Sample  Depicted  as  a  Contour  Map 
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The  common  form  of  the  bivariate  normal   (BVN)   density 
function  is  as  follows: 


f  (X,Y)=• 


2IIQ;CT/(l-p2) 


exp  - 


2(1-PZ) 


'X-yx\*     (X-yx)  (Y-Uy)    /Y-Uy. 


Vr 


a 


y  /  J 


'J 


Since  this  equation  is  that  of  a  probability  density 
function,  its  integral  over  (-oo,  +  ©e>)  must  equal  1.0. 
Consequently,  as  the  standard  deviations  are  increased,  the 
hill  is  flattened.  Additionally,  even  with  small  standard 
deviations,  the  hill  will  remain  less  than  one  unit  high. 
Therefore,  to  be  useful  to  the  modeller,  the  equation  has 
been  modified  (MBVN)  as  follows: 


(1.2)   Z=C  exp  - 


2(l-pz) 


fX-y. 


(X-Uv)  (Y-U„) 


x\      *"  *X'  »-  r-y 


x 


a,ay 


Y-u 


W 


where  C  is  the  peak  elevation  of  the  hill. 

There  are  two  major  advatages  of  the  MBVN  equation. 
First,  it  is  appealing  to  an  observer  as  a  good 
representation  of  many  different  types  of  hills.  Secondly, 
it  has  a  sufficient  number  of  parameters  to  vary  the  shape 
of  a  hill  without  becoming  mathematically  intractable.  The 
parameters  cr^  ,  <7^  ,  and  P  control  the  shape  in  the  XY 
plane,  the  C  parameter  controls  the  height  of  the  hills, 
and///^  and/4y  determine  the  location  of  the  center  of  the  hill 
on  its  map.  (See  Figure  3  for  illustrations  of  the 
parameters) . 
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FIGURES:  Single  Hill  Parameters 
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Ill .   MODEL  METH6D0LGY 

A.   TERRAIN  SIMULATOR 

So  far  in  the  discussion  of  the  new  approach  to 
modelling  terrain,  only  a  single  hill  has  been  discussed. 
This  has  been  important  as  a  mathematical  foundation,  but 
lacking  in  practical  value.  To  be  of  use  in  a  combat 
simulation,  the  hills  must  be  created  collectively  in  a 
size,  quantity,  and  configuration  so  as  to  represent  a 
desired  type  of  terrain.  For  example,  a  user  may  want 
terrain  characterized  by  a  few,  low  rolling  hills;  or 
perhaps  rugged,  mountains  with  peaks  of  widely  varying 
elevations.  Both  can  be  modelled  by  adjusting  parameters. 
The  remainder  of  this  chapter  explains  the  methodology  by 
examining  a  working  terrain  model. 

The  Simulated  Terrain  Model  (SIMTER)  has  two  parts. 
The  first  is  the  main  program  which  creates  terrain  and,  at 
the  discretion  of  the  user,  plots  both  a  three-dimensional 
drawing  and  a  contour  map.  Used  strictly  for  terrain 
generation,  it  can  be  used  as  a  preprocessor  for  other 
combat  models  by  producing  grid  points  and  elevations 
similar  to  those  provided  on  computer  tapes  by  WES.  The 
second  part  of  SIMTER  is  a  movement/line-of-sight  subroutine 
(MOVLOS)  which  moves  a  target  along  specified  routes  across 
the  generated  terrain.  This  part  will  be  discussed  in 
further  detail  in  Chapter  V. 
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In  addition  to  those  input  parameters  used  to  describe 
a  single  hill,  others  are  used  to  aggregate  the  hills  into  a 
map.   These  include  the  following: 

1.  The  dimensions  of  the  battlefield  in  meters. 

2.  The  grid  interval  (e.g.  10  meters  or  100  meters) . 

3.  The  number  of  hills  to  be  created. 

4.  How  much  the  peaks  of  the  hills  are  to  vary. 

5.  How  much  the  spread  of  the  hills  are  to  vary. 

6.  How  many  ridge  lines,  if  any,  are  to  be  created. 

Given  the  input  parameters,  the  model  proceeds  as 
follows. 

1.  The  first  step  of  the  program  is  to  randomly  select  the 
centers  of  mass  for  the  desired  number  of  hills.  This 
operation  is  performed  by  drawing  Uniform  (0,1)  random 
numbers  and  multiplying  them  by  the  size  of  the  battlefield 
in  the  X  and  Y  directions,  thus  producing  an  array  of  paired 
grid  points.  If  ridge  line  are  desired,  then  the  random 
points  in  either  the  X  or  I  direction,  but  not  both,  can  be 
biased  by  drawing  random  normal  deviates  about  a  preselected 
ridgeline  center.  Although  not  a  requirement,  this 
particular  program  makes  a  hasty  plot  of  the  centers  of  mass 
points  for  visual  referrence  of  the  random  process. 

2.  The  second  step  is  to  create  the  dimensions  of  each  MBVN 
hill.  If  the  user  desires  that  all  the  hills  be  of  the  same 
size  and  shape  then  this  step  is  complete.  However,  for 
most  cases  the  height  and  spread  of  each  hill  will  be  varied 
according  to  an  input  parameter  which  is  actually  the 
standard  deviation  for  the  variations.  For  example,  if  a 
user  desires  hills  which  average  100  meters  high,  but  vary 
about  this  average  value  by  approximately  30  meters,  then 
100  becomes  the  mean  value  for  the  peaks  of  the  hills  which 
vary  in  individual  elevations  according  to  normally 
distributed   random   numbers  whose  standard  deviation  is  30. 
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The  same  procedure  is  done  for  the  spread  of  the  hills  in 
both  the  X  and  Y  directions. 

If  only  movement  and  line-of-sight  calculations  are  to 
be  made,  then  MOVLOS  subroutine  is  called  and  the  program  is 
terminated.  A  major  advantage  of  this  approach  to  terrain 
modelling  is  that  no  10  meter,  100  meter,  or  other  grid 
system  is  necessary.  The  terrain  representation  is 
continuous;  therefore,  the  elevation  at  any  point  on  the  map 
can  be  found  without  storing  any  digitized  (discrete)  map 
information.   Only  the  hill  parameters  need  be  stored. 

3.  If  a  matrix  of  grid  values  is  desirable,  as  would  be  the 
case  if  the  program  were  used  in  lieu  of  the  WES  computer 
tapes,  then  the  third  part  of  the  SIMTER  main  program  is  the 
creation  of  an  evenly  spaced  grid  system,  to  include  the 
elevation  at  each  grid  line  intersection.  For  all  the  maps 
and  drawings  in  this  study,  a  100  meter  grid  interval  was 
used. 

Since  the  BVN  distribution  has  some  finite  value 
(elevation)  in  all  directions,  regardless  of  how  far  removed 
from  its  center,  the  elevation  of  each  hill  must  be  computed 
for  each  point  in  the  grid  system.  As  each  hill  is  checked 
for  its  elevation  at  any  selected  point,  it  is  compared  with 
its  predecessor.  If  it  is  lower,  it  is  discarded;  otherwise 
it  is  saved  for  comparison  with  the  next  hill's  elevation  at 
that  same  point.  After  each  hill  has  been  searched,  only 
the  highest  value  is  stored.  It  is  this  value,  along  with 
the  other  maximum  values  at  each  grid  interval,  which  maice 
up  the  terrain  surface.  If  there  is  a  large  number  of 
hills,  the  search  process  can  be  streamlined  by  truncating 
the  MBVN  if  the  height  of  the  density  falls  below  some 
specified  value. 

4.  The  final  part  of  the  program  is  to  provide  the  desired 
output.   This  particular  program  provides  both   printed   and 
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graphical  output.  Both  the  contour  map  and  the 
three-dimensional  drawing  were  produced  by  a  California 
Computers  Compay  CALCOMP  Plotter. 

B.   OSER  INTERFACE 

To  be  of  value  to  the  user,  the  model  must  be  flexible 
and  its  results  realistic.  In  order  to  demonstrate  these 
capabilities  numerous  test  runs  were  conducted.  The  first 
group  of  runs  used  three  different  average  elevations  (100, 
200,   and   300   meters)  and  three  different  numbers  of  hills 

(5,15,25).  The  nine  resulting  terrains  are-  depicted  in 
Appendix  A.  All  maps  and  drawings  represent  a  6,000  meter 
by  3,000  meter  battlefield.  The  second  group  of  runs 
displays   the   results   of   varying  the  seeds  of  the  Uniform 

(0,1)  random  number  generator  which  were  used  to  create   the 
hill  centers   (Appendis   B) .    The  third  group  (Appendix  C) 
demonstrates  the  effects  of   changing   the   spreads   of   the 
hills.   For  all  the  tests  listed  above,  the  contour  interval 
representing  elevation  was  50  meters. 
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IV.   LINE-OF-SIGHT 

A.   THE  PRECISE  APPROACH 

A  major  advantage  of  continuous  terrain  representation 
is  that  line-of-sight  calculations  can  be  made  for  any  pair 
of  grid  coordinates.  This  is  uolike  discrete  terrain  which 
requires  some  approximate  method  such  as  linear 
interpolation.  If  a  simple  "yes"  or  "no"  answer  is  desired 
for  the  question  of  intervisibility  between  any  two  points 
on  the  terrain  surface  then  the  mathematics  consists  of 
solving  the  simultaneous  equations  of  a  straight  line  (0-T 
line)  and  a  curved  surface  in  three  space.  If  a  solution 
exists  then  there  is  not  line-of-sight,  i.e.  the  line 
intersects  the  hills. 


>•  X 


Figure  4:   Line-of-Sight 


The  following  is  a  summary  of  the  procedures 


The  0-T  line  is  parameterized  as  follows: 


(4.1) 

X  =  a  +  b  t 
1           2 

(4.2) 

Y  =  a  +  b  t 
2           2 

(4.3) 

Z  =  a  +  b  t 
3           3 
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The  equations  are  then  substituted  into  the   MBVN 
resulting  in  the  following  expression: 


equation, 


a3+b3t  =  C  exp 


2  (1-p2) 


(4.4) 


(ai+b1t-«x)  (a2+b2-x/y)   /a2+b2t->i/  a 


One  approach   to   solving  this   equation   is   to   take 
natural  logarithms  of  both  sidesr  thus  yielding 


ln(a3+b3t)  =  InC  - 


2 (1-p2) 


[a1+b1t-Ux\ 


\        ax 


(4.5) 


(a1+b]_t-/wx)  (a2  +  b2t-yy) 
-2p + 


'a2  +  b2t-yy 


aa 


V 


Expanding  the  left  hand  side  using  series  expansion  of  a 
logarithm  amd  truncating  after  three  terms  produces  a 
quadratic  expression.  This  can  be  combined  with  the 
quadratic  on  the  right  and  solved  using  the  standard 
quadratic  formula.  Unfortunately,  the  series  has  very  poor 
convergence  with  only  three  terms.  Using  the  average  height 
of  a  hill  as  an  approximate  solution  to  the  series,  less 
than  one  decimal  place  accuracy  can  be  expected.  The 
expansion  used  is  provided  below: 


(4.6) 


(t-a)     (t-a)^ 

In  t  =  In  a  + + 

a        2a 


0<t  6.   2a 


One   approach  to  improving  the  accuracy  of  this  solution  is 
to  solve  the  exact  solution  of  a  similar  0-T  line   which   is 
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horizontal  and  a  hill  which  is  circular  normal,  and  use  the 
results  to  improve  the  convergence  of  the  series  in  just 
three  terms.  An  explanation  of  this  approximation  is 
provided  in  a  later  section. 

Another  way  of  solving  the  problem  of  finding  a 
solution  for  the  parameterized  equations  of  a  line  and  a 
surface  is  to  use  non-linear  programming  techniques.  Using 
the  rapid  convergence  properties  of  a  Newton  search  near  its 
solution,  a  precise  answer  can  be  found  with  minimal 
iterations. 

B.   AN  APPROXIMATE  APPROACH 

If  intervisibility  lines  (0-T  lines)  were  horizontal  and 
the  HBVN  hills  were  circular,  then  the  mathematics  would 
yield  a  quick,  exact  solution.  This  is  done  simlpy  by 
projecting  the  0-T  line  onto  the  XY  plane  and  finding  the 
line  from  the  center  of  the  hill  perpindicular  to  the  0-T 
projection  (See  figure  5)  .  The  intersection  of  these  two 
lines  produces  a  set  of  coordinates  in  the  XY  plane. 
Finding  the  elevation  of  both  the  0-T  line  and  the  hill  at 
this  location  yields  two  numbers  which  can  be  compared  for 
line-of-sight.  If  the  line  is  higher  than  the  hill,  then 
intervisibility  exists. 

Sample  calculations  to  test  the  accuracy  of  this 
approximation  yielded  two  decimal  place  precision  for  0-T 
lines  less  than  15  degrees  from  horizontal  and  MBVN  hills 
whose  spreads  in  the  X  and  Y  directions  did  not  exceed  ten 
per  cent  difference  from  each  other.  Fifteen  degree  angles 
can  be  exceeded  in  two  general  cases:  very  high  hills  and 
very  close  ranges.  In  other  words  the  problem  exists  when 
the  observer  and  target  are  at  significantly  different 
altitudes  when  close  to  each  other.  The  inaccuracy  is 
illustrated  below.  The  curves  represent  two-dimensional 
cross  sections  of  some  part  of  a  hill. 
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Top  View 


Oblique  View 


FIGURES;  Line-of-sight  Calculati 


ON 


22 


LiNE 


EXACT    CASE 


ACTUftL 
INTERSECTION 


REPORTEP 
INTERSECTION 


INACCURATE      CASE. 


Figure  6:   0-T  Line  Intersection  with  Hill  Cross  Section 

An  efficient  solution  to  the  problem  is  to  use  the 
circular  normal  and  horizontal  line  calculations  as  long  as 
the  situation  remains  within  the  confines  listed  above,  and 
use  a  Newton  search  for  those  few  cases  which  fall  outside 
the  desired  accuracy. 
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V .   MOVEMENT/LINE-OF-SIGHT  SaBROaTINE 

A.   MODEL  METODOLOGY 

The  movement/Line-of-sight  subroutine  (MOVLOS)  performs 
three  main  calculations:  movement,  line-of-sight,  and  per 
cent  of  target  area  visible.  To  perform  these  functions 
additional  input  parameters  are  required.  The  first  of 
these  is  the  location  of  a  stationary  observer  or  defender, 
the  second  is  an  array  of  coordinates  which  depicts  the 
route  of  a  moving  target,  and  the  final  input  group  includes 
the  dimensions  of  the  target.  For  this  program  the  target 
is  presumed  to  be  rectangular  so  that  height,  width,  and 
length  are  required. 

The  general  mathematics  of  the  line-of-sight 
calculations  were  explained  in  Chapter  IV.  This  section  is 
designed  to  explain  the  methodology  of  applying  these 
mathematics  to  a  combat  model.  The  first  step  is  to 
determine  the  elevation  of  the  observer  by  searching  for  the 
highest  hill  at  the  observer  coordinates.  This  is 
accomplished  by  substituting  the  observer  XY  coordinates 
into  the  equation  of  each  aBVN  hill.  The  highest  resulting 
functional  value  is  the  elevation  of  the  observer.  This 
procedure  is  likewise  carried  out  for  the  elevation  of  the 
target.  Using  the  three-dimensional  coordinates  of  both  the 
observer  and  the  target,  the  subroutine  forms  the  equations 
for  the  0-T  line  in  three  space.  The  projection  of  the  0-T 
line  onto  the  XY  plane  intersects  at  right  angles  with  lines 
drawn  from  the  centers  of  mass  of  each  hill.  The 
coordinates  of  the  intersections  are  used  to  determine  the 
elevation  of  each  hill  along  the  0-T  line  by  once  again 
substituting  into  the  MBVN  equation  of  the  respective  hills. 
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O-T    LIHE 


Figure  7;   O-T  Line  Projection 

The  elevation  of  the  O-T  line  is  also  computed  at  each 
intersection.  This  is  done  using  the  ratio  of  sides  of 
similar  triangles. 


TARGET 


HTGT 


HUN  = 


<r 


TNTBRSEGTiCN 
DOT  ^ 


DOT 


Figure  8:   Height  of  O-T  Line  at  Intersection 

If  any  hill  is  higher  than  the  O-T   line   along   the   line's 
pathr  then  intervisibility  does  not  exist. 

To   be   of   use   in   a   high   resolution   model,    the 

line-of-sight  subroutine  must  answer  more  than  "yes"  or  "no" 

to  the  question  of  intervisibility.   Since   the   probability 

of  a   hit   (P  )   is   partly   a  function  of  target  size,  the 
h 

program  must  determine  what  effective  area  of  the  target   is 

visible   to   the  firer.   This,  in  turn,  is  a  function  of  the 
target's  dimensions,  the  per  cent  of  area  exposed  above  the 
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terrain,  and  the  angle  it  is  facing  with  respect  to  the 
observer.  Using  basic  trigonmetry,  MOVLOS  computes  the  area 
of  the  target  projected  in  the  direction  of  the  0-T  line. 
Two  possible  target  configurations  are  illustrated  in  Figure 
9.  The  mirror  images  of  these  cases  are  calculated  by  using 
absolute  values  for  the  angles. 


The  dynamics  of  the  program  are  provided  by  the 
movement  portin  of  MOVLOS.  At  discrete  time  intervals,  the 
target  moves  along  the  preselected  routes  toward  its 
objective.  For  all  runs  of  SIMTER  a  one-second  time 
interval  was  used.  At  each  second  an  instantaneous  velocity 
and  a  line-of-sight  are  computed.  If  the  velocity  is  not 
held  constant,  it  is  then  solely  a  function  of  the  slope  of 
the  terrain.  The  slope  is  found  by  taking  the  directional 
derivative  (DEL)  of  the  MBVN  in  the  direction  of  the  target. 
The  directional  derivative  of  f(x,y)  in  the  direction  is 
given  by 


(5.1) 


3f  3f 

DEL  =  —  COS  9  +  --  SIN 
3X  3Y 


For  the  Modified  Bivariate  Normal   the   partial   derivatives 
are-  as_fo Hows: 


3f 
3x 


(1-P2) 


x-u. 


pY-pu 


axay 


expi -- 

2  (1-P2) 


'X-U. 


(5.2) 


(X-«x) (Y-m  ) 
2p 

x  y 


/Y-JU, 


<r 


\~Y     \ 


exp{ 

2  (1-p2) 


IX-v,^ 


(5.3) 


axoy 
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Target  headed  left 
of  observer 


A-B  IS  PROJECTION  OF 
TARGET  SEEN  BY  OBSERVER 


Target  headed  right 
of  observer 


FIGURES:  Effective  Area  of  Target 
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The  output  of  the  Movement/line-of-sight  subroutine 
includes  the  elevation,  coordinates,  speed,  per  cent 
exposed,  and  area  exposed  at  each  one-second  time  increment. 
Additionally,  a  summary  is  provided  at  the  end  of  the 
simulation.   The  information  printed  is  listed  below: 

1.  The  number  and  length  of  intervisiblity  segments. 

2.  Total  distance  traveled  with  intervisibility. 

3.  Total  distance  traveled  without  intervisibility. 

4.  Total  distance  traveled. 

5.  Average  distance  traveled  with  intervisibility. 

6.  Per  cent  of  time  in  which  intervisibility  existed. 
Refer  to  Figures  15  and  16  for  sample  output. 

B.   A  TEST  CASE 

Although  numerous  test  runs  were  conducted,  a  single 
random  sample  was  initially  chosen  so  that  an  evaluation 
could  be  made  of  the  model's  usefulness  in  a  real  situation. 
The  intent  was  to  conduct  the  test  in  a  manner  similar  to 
the  actual  projected  use  of  the  model.  In  this  regard,  the 
following  steps  occurred: 

1.  The  military  analyst  selected  a  "type"  of  terrain. 

2.  The  user  applied  the  terrain  description  to  the  model  by 
adjusting  the  input  parameters  (Figure  10) . 

3.  The  user  executed  SIHTER  without  the  HOVLOS  subroutine 
in  order  to  obtain  a  graphical  representation  of  a  single 
random  realization  of  terrain  (Figures  11  and  12). 

4.  The  analyst  selected  a  location  for  the  defender 
(observer) ,  and  routes  for  the  attacker.  The  routes  were 
chosen  based  on  military  judgement.  For  this  case  routes 
were  based  on  maximizing  cover  while  attempting  to  maintain 
good  visual  surveillance  of  the  battlefield  (Figures  13  and 
14)  . 

5.  The  analyst  provided  maximum  and  minimum  speeds  for  the 
attacker.  For  this  run  the  minimum  speed  was  two  meters  per 
second  and  the  maximum  speed   was   ten   meters   per   second. 
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These  speeds  occurred  at  a  45  degree  positive  slope  and  a  45 
degree  negative  slope,  respectively.  Since  the  relationship 
was  linear,  level  terrain  velocity  was  six  meters  per  secon. 
6.  The  user  resubmitted  the  program,  this  time  without 
creating  graphics  but  including  MOVLOS  instead.  The 
resulting  output  (Figures  15  and  16)  compared  logically  with 
what  appeared  on  the  terrain  maps. 
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Figure  11:   A  Sample  Terrain 
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Figure  14:   Sample  Terrain  with  Routes 
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C.   VERIFICATION 

To  ensure  that  the  model  was  responding  to  all  changes 
in  parameters,  a  series  of  verification  runs  were  conducted. 
It  was  not  the  intent  of  these  replications  to  establish  the 
sensitivity  of  line-of-sight  to  variations  in  terrain.  The 
latter  was  the  subject  of  the  previously  mentioned  VRI 
study. 

The  first  series  of  runs  consisted  of  varying  the 
number  of  hills  (5,25,  and  25)  and  the  height  of  the  hills 
(100, 200, and  300).  Since  the  same  random  number  generator 
seeds  were  used  for  all  nine  runs,  the  15  hills  included  the 
original  5,  and  the  25  hills  included  the  original  15.  This 
made  it  possible  to  use  the  same  routes  for  all  replications 
without  violating  tactical  doctrine  (Appendix  D) .  A  matrix 
of  the  results  is  presented  in  Table  1 . 

A  second  group  of  tests  was  conducted  by  varying  the 
maximum  and  minimum  speeds  of  the  attacker  over  the  same 
range  of  slopes  (-45, +45).  The  initial  (2,10)  meters  per 
second  was  run  against  (4,16)  meters  per  second.  The 
results  yielded  no  significant  difference  between  the  two 
runs  except  that  the  faster  target  was  intervisible  95  out 
of  501  seconds  while  the  slower  was  intervisible  160  out  of 
836  seconds  (both  1 955)  . 


37 


2 


on 
on 


o 
o 

I— • 

I— 


a: 


a> 

f— 

.a 

•^ 

O) 

to 

<u 

u 

•^ 

O) 

c 

> 

00 

O 

O 

CM 

^~ 

r— 

CO 

CO 

CO 

fO 

<T3 

4_ 

rv. 

r>» 

CO 

0 

«3- 

r^. 

CO 

«3- 

s_ 

+-> 

Oi 

'O 

co 

r— 

CO 

r>. 

•=*- 

CO 

«3- 

CM 

dj 

CO 

4-> 

CM 

CO 

CM 

> 

•^ 

c 

<  0 

t— 1 

ai 

r— 

-Q 

«4- 

•r- 

O 

CO 

CO 

4-> 

S_ 

> 

c 

CM 

r— 

CO 

CO 

CO 

co 

,— _ 

,_ 

CO 

CD 

s- 

a» 

0 

-Q 

O) 

E 

0 

E 

4-> 

C7> 

CM 

3 
OJ 

c 

1—1 

a> 

CO 

a  -0 

c 

a> 

«3 

s_ 

CVJ 

O 

r— 

r^. 

10 

CM 

co 

CO 

CO 

■(-» 

O) 

• 

• 

• 

• 

• 

• 

• 

CO 

> 

P>. 

O 

CNJ 

CO 

1 — 

O 

CO 

00* 

0 

•r*" 

0 

00 

O 

cn 

CO 

*r 

CM 

O 

0 

■"" 

<3« 

rD 

<u 

+-> 

3 
0 

OS 

<c 

CQ 

0 

«=c 

CO 

O 

< 

0Q 

c_> 

<u 

l^~ 

-Q 

•r™ 

(D 

CO 

0 

U 

•F— 

CD 

c 

> 

0 

O 

■ — 

CO 

r— 

CO 

O 

CO 

CT> 

03 

ra 

s_ 

co 

r*» 

r- 

CO 

r>» 

CO 

r-» 

CO 

S_ 

4-> 

O) 

c£> 

on 

CSJ 

(■>•» 

^r 

fssm 

co 

«tf- 

CM 

<D 

00 

+J 

CVJ 

10 

CO 

co 

fmm 

> 

•r- 

c 

<  Q 

t— 1 

<u 

f— 

-O 

4- 

•^- 

0 

00 

co 
+■> 

S- 

> 

c 

CSJ 

n- 

C\J 

CO 

,— 

CM 

,_ 

,__ 

CO 

a> 

i_ 

OJ 

0 

-Q 

O) 

E 

0 

E 

a) 

c 

1—1 

CD 

a> 

CO 

u  -0 

c 

aj 

03 

s- 

ir> 

O 

CO 

00 

00 

CO 

cr» 

CM 

CO 

4-> 

aj 

CO 

> 

CO 

O 

CSJ 

r^ 

<s- 

CM 

co 

CO 

,_ 

•r- 

0 

00 

O 

<y> 

CO 

CO 

CO 

,— _ 

Q 

0 

*"~ 

^  Z3 

<U 

4-> 

3 
O 

<C 

ca 

<_> 

<c 

CQ 

C_> 

<C 

CO 

c_> 

to 

a, 

4-4 

o 

a 

o 

•H 
+J 

ccS 
> 

i—( 

w 

03 

CO 


4-4 

o 

E 
3 


00 

fi 

•H 
X 

u 

> 

4-4 
O 

CO 


3 

CO 

CD 


1— I 
03 


co 


CO 


CO 
CM 


S3cjniV3J  NIVHcd31  JO  cd33WnN 


38 


A  third  series  of  runs  consisted  of  varying  the  spread 
of  the  hills  in  the  initial  test  case,  while  holding  all 
other  parameters  constant.  As  was  expected,  intervisibility 
decreased  as  the  hills  were  widened. 


SPREAD 

PEK  CENT  INTEEVISIBLE 

300 

27.0 

400 

19.4 

500 

7.7 

Table  2:   Effect  of  Varying  Spread  on  Intervisibilty 

Other  runs,  such  as  varying  the  location  (seeds)  of  the 
hills  while  holding  routes  constant,  were  executed  but  not 
included  in  this  report. 
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71. 


If  exact  terrain  modelling  is  not  required  in  a  combat 
simulation,  then  representative  terrain  can  be  created  using 
a  modified  bivariate  normal  distribution.  Since  there  is  no 
requirement  for  survey  or  photographic  interpretation  in 
order  to  mathematically  model  terrain,  this  approach  is 
significantly  less  costly  and  time  consuming  than  digitizing 
terrain.  An  additional  advantage  of  mathematical 
representation  is  that  replications  of  a  "type"  of  terrain 
can  be  randomized,  thereby  improving  the  statistical  level 
of  confidence  in  a  combat  model  output. 

The  parameterized  terrain  is  continuous;  therefore,  the 
elevation  at  any  location  is  exact  and  not  a  linear 
interpolation  between  discrete  points.  This  makes  it 
possible  for  line-of-sight  calculations  to  be  more  accurate. 
On  the  other  hand  if  the  program  is  to  be  used  as  a  terrain 
preprocessor  for  a  high  resolution  combat  model,  then  SIMTEB 
can  produce  a  digitized  output  from  the  parameterized 
representation. 

The  results  of  the  tests  indicate  that  the  parametric 
representation  of  terrain  is  both  useful  and  realistic.  In 
view  of  these  conclusions  the  following  action  is 
recommended: 

1.  Parameterized  terrain  should  be  run  against  digitized 
actual  terrain  in  a  high  resolution  model  such  as  DYNTACS  or 
AMC  74.  Sufficient  replications  of  the  randomized 
parametric  terrain  should  be  conducted  in  order  to  establish 
steady  state  results. 

2.  Distributions  other  than  the  bivariate  normal  should  be 
examined,  e.g.  the  beta. 

3.  The  feasibility  of  representing  actual  terrain  with  the 
MBVN  should  be  examined  by  preselecting  centers  of  hills  as 
they  appear  on  a  map. 


40 


4.  The  SIMTER  simulation  should  be  evaluated  for  use  as  a 
mobility  model.  Houtes  can  be  selected  or  readily  changed 
commensurate  with  vehicle  performance. 
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APPENDIX  A 
MATRIX  OF  INPUT  VARIABLES 

The  graphs  included  in  this  appendix  illustrate  variations 
in  number  of  terrain  features  (N)  and  average  height  of  hills 
(PEAK).   The  graphs  are  presented  in  pairs;  the  first  being  the 
contour  map  and  the  second  being  its  three  dimensional  repre- 
sentation.  The  order  of  variation  is  listed  below. 
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APPENDIX  B 
VARIABLE  RANDOM  NUMBER  SEEDS 

The  graphs  included  in  this  appendix  have  the  same  constant 
"nput  parameters.   The  only  variable  is  the  seed  for  the  random 
lumber  generator.   As  with  the  previous  appendix,  the  graphs 
appear  in  pairs;   the  contour  map  first,  followed  by  its  three 
dimensional  representation. 
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APPENDIX  C 


VARIABLE  SPREAD  OF  HILLS 


The  graphs  presented  in  Appendix  C  have  constant  input 

parameters  except  for  the  spread  of  the  hills  (XMEL,YMEL). 
The  three  variations  are  100,200,  and  300  meters,  respectively 
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APPENDIX  D 
ROUTE  SELECTION 

The  contour  maps  in  Appendix  D  are  the  same  as  those 
in  Appendix  A  except  for  the  overlayed  routes.   These  graphs 
were  used  to  develop  input  data  for  the  1 i ne-of-s ight 
subroutine. 
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